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Abstract. We calculate a pump current in a classical two-state stochastic chemical 
kinetics by means of the non-adiabatic geometrical phase interpretation. The two-state 
system is attached to two particle reservoirs, and under a periodic perturbation of the 
kinetic rates, it gives rise to a pump current between the two-state system and the 
absorbing states. In order to calculate the pump current, the Floquet theory for the 
non-adiabatic geometrical phase is extended from a Hermitian case to a non-Hermitian 
case. The dependence of the pump current on the frequency of the perturbative kinetic 
rates is explicitly derived, and a stochastic resonance-like behavior is obtained. 
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1. Introduction 

Some of classical stochastic processes have been investigated via analogies with quantum 
mechanics, and such analogies give us useful analytical methods for investigating various 
physical quantities and phenomena in classical stochastic processes. For example, 
asymmetric simple exclusion processes have been studied a lot by means of the free 
fermionic representation or the quantum spin formalism [1]. For the reaction-diffusion 
processes, it has been shown that the second-quantized representation and field theoretic 
approach are useful [2-5]; deeper understandings of the non-equilibrium behavior has 
been achieved with the aid of the renormalization group method [6]. Recently, the 
second-quantization method has been applied to stochastic switching phenomena in 
gene regulatory networks, with the aid of a variational method [7-10]. 

Recently, Sinitsyn and Nemenman [11, 12] have been studied a pump current in 
stochastic chemical kinetics. They studied a classical two-state stochastic system in a 
sea of substrates and products (absorbing states); the system can be interpreted as a 
single Michaelis-Menten catalyzing enzyme or as a channel on a cell surface. Under 
a periodic perturbation of the kinetic rates, the classical system gives rise to a pump 
current between the absorbing states. By means of the expression for the full counting 
statistics of transitions among the absorbing states, they have clarified that there is a 
relationship between the Berry phase and the pump current, i.e., the pump current can 
be separated into two parts: one of them is a classical current, and the other is indeed 
a geometrical one which depends only on the contour in the parameter space. Sinitsyn 
and Nemenman have calculated the geometrical phase in the adiabatic case, in which 
the oscillations of the perturbative kinetic rates are very slow [11]. It has been clarified 
that the pump current is linearly proportional to the frequency in the low frequency 
regime. 

In contrast with the low frequency case, such a pump current is expected to show 
a power-law decay behavior with respect to the frequency in the high frequency regime 
[13-15]. In addition, there should be a single maximum in the pump current, which is 
a similar behavior to the stochastic resonance. While this behavior has been obtained 
by another approach for a system slightly different to the Sinitsyn and Nemenman 
one [14], it is expected to be possible to obtain the behavior of the pump current by 
the 'geometrical interpretation.' Furthermore, such studies for classical systems from 
the viewpoint of the quantum mechanical formulation might become a key step towards 
deeper understanding of the classical stochastic systems. 

In the present paper, we analyze the pump current in the 'non-adiabatic' case. 
In order to evaluate the pump current in the non-adiabatic case, we will extend the 
usual Floquet theory, which has been discussed by Moore [16], to a non-Hermitian 
case. In addition, the Floquet Hamiltonian cannot be block diagonalized in our case, 
and hence a perturbation theory is applied. From these calculations, we can derive an 
explicit dependence of the pump current on the frequency of the perturbative kinetic 
rates, which shows a linear dependence on the frequency in the low frequency regime 
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(adiabatic case), and a power-law decay in the high frequency regime. 

The construction of the present paper is as follows. In section 2, we explain a 
model for the pump current, which was introduced by Sinitsyn and Nemenman [11]. 
The general formalism for the calculation of the pump current is given in section 3. 
Section 4 is devoted to an extension of the Floquet theory developed by Moore [16] to a 
non-Hermitian case. In section 5, we exhibit an explicit calculation of the pump current 
with the aid of perturbation theory. Section 6 gives concluding remarks. 



2. The model 



Sinitsyn and Nemenman [11] have introduced a candidate minimal models for the pump 
current. The classical chemical kinetic system is constructed as follows: 
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The system consists of three parts. The container can contain either zero or one particle 
in it. When the container is filled with one particle, the particle can escape from the 
container by jumping into one of the two absorbing states: the left reservoir or the right 
one. In contrast, when the container is empty, either of the absorbing states can emit a 
new particle into the container. Here, we define the pump current as the mean particle 
current J from the container into the right reservoir. 

For simplicity, in the present paper, we consider the following kinetic rates: 

£;_i = k 2 = a, 

k x = b + Rcos(u>t), (2) 
A;_2 = b + Rsm(u)t), 

where a and b are positive real numbers, and R is the amplitude of the perturbative 
oscillation. The above kinetic rates mean that only two kinetic rates (hi and fc_ 2 ) 
oscillate with time at a frequency u. For the above kinetic rates, there is no classical 
L — > R current, as discussed in reference [11]. Hence, if there is a current, it is called a 
pump current, which stems from the geometrical term. 



3. Geometrical phase and pump current 



We here give a brief summary of the discussion given by Sinitsyn and Nemenman [11], 
and extend their formulation to one useful for the non-adiabatic case. 

Let Pf and P c be the probabilities for filled and empty container respectively. When 
we define a state of the system by 



P(*) 



the master equation for the time evolution of the system is written as 



df 



P(*) = 



fci + /c_2 -k-i - k 2 
—ki - k_ 2 k-i + k 2 



Pit). 



(3) 



(4) 
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Note that P e = 1 - P f . 

In order to investigate the pump current between the container and the right 
reservoir, we define P n as the probability of having n net transitions from the container 
into the right reservoir during time t. Using the probability P n , the probability 
generating function for the number of transitions is written by 

oo 

Z( X ) = e S(x) = £ p n= s e isx , (5) 

s=— oo 

where x is called the counting field, and S(x) the full counting statistics. Formally, the 
generating function Z(x) is given by [11] 

Z( X ) = l1f(e-&*< x * ) *)p{t ), (6) 



H( X ,t) 



(7) 



where 1 is the unit vector, T stands for the time-ordering operator, and 

kx + /c_2 —k-i - k 2 e lx 
-k\ - fc_ 2 e~ lx k-i + k 2 

The derivatives of S(x) gi ye cumulants of P n , e.g., (n) = —idS(x)/dx\ x =o- We here 
denote the period of the rate oscillation by T(= 2tv/u). When we assume that a state 
p(t) gives a cyclic evolution of the system, i.e., p(T) = e K ^p(0), the pump current is 
given by J pump = (—i/T)dK(x)/dx\ x =Oi because of the normalization condition of the 
probability (l t p(0) = 1). 

In order to calculate the full counting statistics (or K(x)), we focus on the fact that 
the probability generating function is related to the solution of the following differential 
equation: 

± p (t) = -H( X ,t)p(t). (8) 

When we replace the time-evolution operator H(x,t) by H = — iH(x, t), we obtain the 
following Shrodinger-like equation: 

i^0(t) = Hftt). (9) 

Hence, the problem of the calculation of the full counting statistics (or K(x)) is replaced 
by the problem of the evaluation of the geometrical phase in equation ([9]). We here 
note that the 'Hamiltonian' if is a non-Hermitian operator, so that we need a slight 
modification of the usual calculation for the geometrical phase. 



4. Extension of Moore's discussion to a non-Hermitian case 

In order to perform the calculation in the non-adiabatic cases, one can use the Floquet 
theory discussed by Moore and Stedman [17]. Choutri et al. [19] have extended the 
discussion given by Moore and Stedman to non-Hermitian cases, and calculated the non- 
adiabatic geometrical phase for a special case [19]. However, in the above formulations it 
is difficult to give an explicit expression for the geometrical phase for general cases. On 
the other hand, Moore [16] have combined the above Floquet theory with the discussion 



The stochastic pump current and the non-adiabatic geometrical phase 



5 



given by Shirley [18], and obtained a useful formulation for the non-adiabatic geometrical 
phase in the Hermitian case. We will give the extension of the Moore's discussion to 
non-Hermite cases, and hence the following discussion is mainly based on the Moore's 
one [16]. 



4-1. The Floquet theorem 

At first, we give a brief summary for the Floquet theorem [17,20,21]. A non-singular 
matrix A is called a fundamental matrix of the system of ordinary differential equations 
x{t) = W{t)x{t) if A{t) = W{t)A{t), where x{t) and A(t) mean the time derivative 
of x(t) and A(t), respectively. Let W(t) be periodic in t. Then, the Floquet theorem 
ensures that there exists a B(t) periodic in i and a constant C such that A(t) = B(t)e ct . 



4-2. The non-adiabatic geometrical phase and the Floquet theory 

For a periodic non- Hermitian Hamiltonian H with period T, i.e., H(t + T) = H(t), 
the non-unitary time-evolution operator U(t) is decomposed into the Floquet product 
form as U(t) = V(t)e lMt . Due to the Floquet theorem, the non-unitary operator V(t) 
has a period T, and M is a time-independent operator. In addition, U (t) is the unique 
fundamental matrix satisfying U(0) = I because U(t) is the time-evolution operator. 
In contrast to the Hermite Hamiltonian case discussed in reference [16], the T-periodic 
matrix V(t) is non-unitary, and M is non-Hermitian. Because of the non- Hermitian 
property of the Hamiltonian, we here introduce the non-unitary time-evolution operator 
of the periodic 'adjoint' Hamiltonian H^(t) [19]. The operator is written as U(t), and 
we decompose U(t), in a similar way to U(t), as U{t) = V(t)e lMt . Choutri et al. [19] 
have shown that the non-adiabatic geometrical phase is given by 

ia = f T (Mo)\^v\Mo))dt, (io) 

JO 

where |0«(O)) are the right eigenvectors of U(T), and \4> a (Q)) are those of U(T). Each 
\4> a (0)) is also an eigenvector of M, because U{T) = e lMT . We here note that (0 Q (O)| is 
not the conjugate of |0 Q (O)), but that of |0 Q (O)). Finally, the pump current is evaluated 
from 

T .... ld_ _ u_d_ nn 

■'pump rr o la a la- 

T dx 2tt dx 

However, it is difficult to obtain V(t) and M directly in general; some additional 
analytical treatments are needed in order to calculate the geometrical phase j a . 



4-3. Geometrical phase in terms of the Shirley's fundamental matrix 

There are two fundamental matrices of particular importance for calculating the 
geometrical phase [16]. The first is the time-evolution matrix U(t) which is discussed 
above. The second is the fundamental matrix F(t) = P(t)e lQt , which is discussed 
by Shirley [18]. The fundamental matrix F(t) has suitable properties for the Fourier 
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analysis, which enables us to calculate V(t) and M explicitly. In the Hermitian case 
discussed in references [16,18], the matrix Q is real and diagonal in some convenient 
basis. In the non-Hermitian case, the matrix Q is not diagonalized in general, and hence 
all discussions would be modified by using the Jordan canonical form (but we do not 
use the Jordan canonical form explicitly because the matrix Q can be diagonalized in 
our case, as shown later.) 

Henceforth, we will sometimes omit the time dependence of operators (matrices) 
for notational simplicity (such as U = U(t)). Since U and F are both fundamental 
matrices of H, there exists a constant invertible matrix X with U = FX [16]. Because 
U(0) = I, we have 17(0) = F(0)X = I. Hence, we obtain X = F(O)" 1 and 

U(t) = F(t)F(0)^ x = P(t)e iQ< P(0)- 1 = P(t)P(0)^ 1 e iP(0)QP(0)_lt . (12) 

We can here make the identifications V(t) = P(t)P(0)" 1 and M = P(0)QP{0)~\ When 
we set the eigenvectors of Q as \a), the eigenvectors of M, |0 Q (O)), are related to |a) 
by |0 Q (O)) = P(0)|o;). We perform the same calculation for the time-evolution operator 
U(t), and obtain V(t) = P^PiO)- 1 and (0 a (O)| = (5|P(0) t . The non-adiabatic 
geometrical phase is therefore rewritten as 

la = ^(Mm [P(O)^ 1 p(t)tp(t)p( O )^|0a(O))dt 

= [ T (a\iP(tyP(t)\a)dt. (13) 

JO 

4-4- The Floquet Hamiltonian 

So far, we have shown that it is possible to write the non-adiabatic geometrical phase in 
terms of the fundamental matrix F(t) = P(t)e 1( ^* even in the non-Hermitian case. The 
next step is the calculation of the explicit form of P(t). For that purpose, Shirley has 
introduced the Floquet Hamiltonian Hp as follows. We introduce the Fourier expansion 
of the Hamiltonian H as 

oo 

H= Yl H (n) e inuJt , (14) 

n=— oo 

where the Fourier components are defined by 

i7 (n) = (1/T) [ T He- inuJt dt. (15) 
Jo 

The Floquet Hamiltonian is defined by the matrix elements [18] 

(a', n\ H F m) = H^ m) + nu8 a >p5 nm , (16) 

where \a,n) is an orthonormal basis for the matrix representation of H-p. The index a 
represents a spatial part and the index n merely represents a Fourier components. In 
our case, we have only two spatial parts (as we will see it later), and denote them by 
+ or — , i.e., a, f3 e {+,—}• The \a,n) is called a 'Floquet state' [18]. Note that H is 
non-Hermitian, so (a',n\ is not a conjugate state of \a,n), but the corresponding left 
orthonormal basis; we added the prime on a in order to clarify this fact. 
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4-5. The relationship between the Floquet Hamiltonian and the Shirley fundamental 
matrix 

Like for the case with the Hamiltonian H, we define the Fourier expansion of P(t) as 

00 

p= P (n) e inu ', (17) 

n=— 00 

Shirley [18] has shown that the Fourier components of P have a simple relationship 
between the Floquet Hamiltonian in the Hermitian case. Let Hp have eigenvectors 
|£a,n) an d eigenvalues £ a ,n- In the Hermitian case, the matrix elements of Q and the 
Fourier coefficients of P are given by Q a p = —s a ^5 a p and P^fl = (a',n\ep t o) [18]. 

In the non-Hermitian case, the matrix Q is not diagonalized in general, so that 
we cannot continue the same discussion as for the Hermitian case. However, when 
the eigenvectors of the Floquet Hamiltonian, \e a>n ), are not degenerate, the matrix 
Q is diagonalized and we obtain the same expression for Q a p and Pf$ as for the 
Hermitian case. Actually, in our case, the eigenvectors of the Floquet Hamiltonian 
is not degenerate, as shown in section 5 [22]. 

In order to derive the Fourier expression of P\ we introduce another definition of 
the Fourier transformation 

£f(») = (l/T) [ T He inuJt dt, (18) 



and expend P^ as 

00 

pt= pW e - iwt . (19) 

71= — CO 

By using the definition of the Fourier components in equation ffl8|) . other Floquet 
Hamiltonian and Floquet states are introduced: 

(a7n| H F j3', m) = H { ^7 m) - nuj5 a ^5 nm . (20) 

The matrix elements of [P^) a /3' are given by (£ a> o\(3',n), where (£ a ,o| are the 'left' 
eigenvectors of the Floquet Hamiltonian Hp. (Unlike for the case of Hp, we add 
the prime on a for the right eigenvectors of the Floquet Hamiltonian Hp.) From the 
orthogonality relations J T exp[i(n — m)ujt)dt = T5 nm , we obtain 

00 

Lptp= (-nw)P (m)t P (n) e i(n - mM , (21) 

n,m=— 00 

and finally, the geometrical phase is expressed by means of the Fourier components of 
P and P^ as follows: 



la = { « 



]T (-27in)P^P^ n 



a . (22) 
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Table 1. A part of the Floquet Hamiltonian Hp. 
Hamiltonian is block diagonalized. 



When R = 0, the Floquet 



\+,n + l) 



— ,n - 



1) 



-,n> 



I- n) 



|+,n-l> 



|- n-l> 



-2i6 + (n + 
6(i + ie-^) 

He" 1 *) 





a(i + ie' x ) 
-2ia + (n + l)w 







±(-iR-R) 
|(iii + ile-^) 
— 2i6 + no; 
b(i + ie" 1 *) 
±(-iR + R) 






a(i + ic lx ) 
— 2ia + naj 










i(-LR-R) 
^(iR + Re-^) 
-2ib + (n - l)w 

6(i + ie- ix ) 








a(i + ie' x ) 
-2ia + (n — l)a; 



5. Perturbation calculation for pump flux 

5.1. Perturbation for the Floquet eigenvectors 

In order to calculate the Fourier components of P and P\ we need to calculate the 
explicit expression of the Floquet Hamiltonian and its eigenvectors. While it is easy 
to calculate such eigenvectors in the block diagonalized Floquet Hamiltonian as in 
the case in [19], the Floquet Hamiltonian in our case cannot be block diagonalized. 
Table [T] gives a part of the Floquet Hamiltonian Hp. When R = 0, it is easy to see 
that the Floquet Hamiltonian is block diagonalized. Hence, we here use the standard 
perturbation techniques, though a little modification is needed to calculate the perturbed 
states for the left eigenvectors of the Floquet Hamiltonian Hp. 

When R = 0, we can show that the infinite-dimensional Floquet Hamiltonian is 
block diagonalized with typical block 

— 2ib + nuj ia + iae' x 
ib + ibe~ lx — 2ia + nuo 

in the basis {|+, n), | — , n)}. We define the perturbative part H F as the rest of the 
original Floquet Hamiltonian with R ^ 0. 

From the Hamiltonian ( l23i) . we obtain the eigenvalues 

ei 0) n = -ib -ia± ie~ lx ^e ix {e ix b + a)(b + e ix a) + nu, (24) 
and the corresponding right eigenvectors 

|ei 0) n ) = z[-e ix b + e ix a ± ^e ix (e ix b + a)(b + e ix a)] \+, n) + z(l + e ix )b\-, n) . (25) 
The left eigenvectors of H Q are obtained from 

= z'[-e~ ix b + e~ ix a ± Je~ ix {e- ix b + a){b + e~ ix a)} (+', n\ 



(23) 



+z'{l +e- ix )a(-',n\ . (26) 

z and z' are normalization constants satisfying (e^i n \ £ ^ln) = ^a'p^nm- 

From the usual perturbation theory, we obtain the first order correction for the 
state |^^_^^ as 

(JO) | Rl | (0) » , (0) i rr, | (0) v 

u (l) v _ \ E +',n+l\ rl F\ E +,n/ I (0) v . \ E -',n+l \ rl F \ E +,n/ , (0) v 
-,n) — _( ) (0) \ £ +,n+l) + (o) (0) \ £ -,n+l) 



£ +,n £ +,n+l £ +,n £ — ,n+l 
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\ t +',n-l\ rl ¥\ t +,nl I (0) v . \ t -',n-l\ rl F\ t +,n/ , (0) v /n-x 
(01 \ £ +,n-l) + f01 fOl l £ -,n-l/- I/' J 



(0) v . \ c - , ,ra-l|-'-'Fl c +,ra/ | (0) 
if I .(0) ^+,n~V + (0) _ (0) 

+,n— 1 fc +,n fc — n— 1 



The first order correction for |e_ n ) is calculated in a similar manner. The matrix 
»(«) 

= <«>l (l4°o) + l4'o>) • (28) 



elements of P^l are therefore given by 



For the Floquet Hamiltonian Hp, we perform the similar calculation, but in the 
present case, we need the first order correction of the left eigenvectors (e+J. The 
Floquet Hamiltonian Hp has different Floquet states to Hp, and the diagonal part of 
the matrix representation of Hp is given by 

(+,n\Hp\+',n) = —2ib — nu, 



(29) 

(— , n\Hp\— ', n) = — 2ia — nuj. 

The other matrix elements of the Floquet Hamiltonian Hp are the same as those of the 
Floquet Hamiltonian Hp. Hence, we need a slight modification for the case of Hp as 
follows: 

—2ib — nu ia + iae lx 
ib + i&e~ lx — 2ia — nu 



Ha 



(30) 



4°l = ~ib - ia ± ie" ix v /e ix (e ix 6 + a){b + e ix a) - no;, (31) 

|£tf? n ) = ^[-e ix 6 + e ix a ± ^(e^fo + a) (6 + e*a)] |-PTn) + z(l + e lx )b\^n), 

(4 0) n| = z'[-e-' lx b + e'^a ± ye^e"^ + a) (6 + e" ix a)](+^,| + + e - ix )a(^n]. 

The first order correction for the (e+ n | is calculated from 

/jr(0) |7f/ |~{0) v ,~(0) ,77, ,~<0) v 

/~(1) I N^ + .nl-^Fl^+'.n+l/ ,~(0) i \ fc + ,n |-"F F -',n+l/ ,~(0) 
\ £ +,n| — _( ) _ ~{0) \ £ +,n+ll H ~<0) _ ~(0) \ £ -,n+ll 

£ +,n £ +,n+l £ +,n £ -,n+l 

/=(0) iTf/ i~(0) v ,~{0) ,7f, - J0) v 

<F+,ra|- n FF+',n-l/ /~(0) I \ £ +,n|-"FF-',n-l/ /~(0) i /qq\ 
+ ~(0) _ ~{0) \ £ +,n-l \ + ^o) _ ~(0) \ £ -,n-ll> l^J 

£ +,n £ +,n— 1 £ +,rt £ — ,n— 1 

and the matrix elements of J>( n )t are given by 

[pWt U = ( ( ^o)| + <g<J) [) |/3^n>. (34) 

5.i?. Explicit expressions for the pump current and some discussions 

From equations (122]) . (!28|) and (I34l) . we finally obtain the geometrical phase 7+ and 7_ 
after tedious calculation. We here note that there are two cyclic states corresponding to 
the spatial part + and — . While there are two cyclic states, we expect that the pump 
flux to be given by 7+. The reason is as follows. Since the coefficients of the states 
|+, n) and |— ,n) are expected to be related to the probability p(t), these coefficients 
should be positive when x = 0. From the definition of the kinetic rates, we have 
positive parameters a and b (a,b > 0), so that we expect that the eigenvector with 
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the + spatial part gives a suitable geometrical phase for the classical pump current. 
(For the eigenvector with the — spatial part, one of the coefficients of |+, n) and |— , n) 
in equation ( 125]) becomes negative when x = 0.) When 7+ is adopted, we obtain the 
following pump current using equation (TTTT) : 

Jpump = l~ a R \{a + bf + ^ (35) 
and in particular, for the adiabatic case (uj <C 1), 

.(adiabatic) b 2 

Vmp " i 6a(a + 6) 2^ w - W 

We note that the geometrical phase in equation fl22|) consists of the two eigenvectors 
I a) and (5|, and hence the first-order correction in R of the Floquet eigenvectors gives 
the second order in R of equations (1351) and (|36|) . This is also consistent with the result 
in [11] in the small R regime [23]. In addition, the dependence of the pump current on 
the frequency uj is consistent with the pump current in [11]. Hence, we consider that 
the eigenvector with the + spatial part gives a suitable pump current in our formalism. 

For the case with uj 1, equation (1351) shows the power- law decay ~ l/o;, which 
is expected by the other discussions for the pump current [13-15]. Furthermore, we 
have a stochastic-resonance-like behavior with respect to the frequency uj: there is a 
single maximum of the pump current at uji = 2(a + b). This behavior is the same as 
those obtained experimentally and theoretically [13-15]. We therefore conclude that the 
discussion based on the geometrical phase for the pump current gives the pump current 
adequately. 



6. Concluding remarks 

In the present paper, we calculate the pump current of a classical chemical system 
by means of the geometrical phase. First, we extend the Floquet theory to the non- 
Hermitian case, which is applicable for the Floquet Hamiltonian with non-degenerate 
eigenstates [22]. Secondly, we obtain an explicit expression for the pump current with the 
aid of a perturbation technique. While we analysed only a specific case, our formulation 
is applicable for other examples [24]. One of the advantages compared to the known 
techniques is that our formulation gives a formal recipe for the calculation of the pump 
current. Combining the Floquet theory and the perturbation techniques, the explicit 
dependence of the pump current on the frequency uj is derived. In addition, it is 
possible to calculate the fluctuations of the pumping effects (higher cumulants) using 
the derivatives of the geometrical phase [11]. 

The derived pump current is consistent with the previous results which have been 
obtained by the other approaches: the pump current depends on the frequency of the 
perturbative kinetic rate linearly at low frequency regime. In addition, there is a power- 
law decay of the pump current with respect to the frequency in the high frequency 
regime. Note that we use the perturbation method in our calculation, so that the final 
result is provided in the limit of weak perturbations. In order to get quantitatively 
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good results even in the large R cases, it would be needed to consider higher orders of 
R or to use other analytical treatments. However, the dependence of the pump current 
on the frequency u is consistent with the previous results, and we believe that our 
formulation would be useful for further discussions about the geometrical phase and the 
pump current. 

These treatments based on the geometrical phase would give a unifying theory 
for such pumping effects. In addition, the analogy with such geometrical phase in the 
quantum mechanics might give us a deeper understanding of the classical system. We 
expect that the physical interpretation of the eigenvectors of the Floquet Hamiltonian 
might give us more information about the current, and it would be possible to obtain 
deeper understanding of the stochastic resonance-like behavior by considering it with 
the aid of the geometrical interpretation. 
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